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Abstract. A geometric Brownian motion with delay is the solution of a stochastic 
differential equation where the drift and diffusion coefficient depend linearly on the 
past of the solution, i.e. a linear stochastic functional differential equation. In this 
work the asymptotic behavior in mean square of a geometric Brownian motion with 
delay is completely characterized by a sufficient and necessary condition in terms of 
the drift and diffusion coefficients. 



1. Introduction 

Geometric Brownian motion is one of the stochastic processes most often used in 
applications, not least in financial mathematics for modelling the dynamics of security 
prices. More recently, however, modelling the price process by a geometric Brownian 
motion has been criticised because the past of the volatility is not taken into account. 

The geometric Brownian motion is the strong solution of the stochastic differential 
equation 

dX(t) = bX(t)dt + aX{t)dW{t) fort^O, 

where b and a are some real constants. If we wish that the dynamics of the process 
X at time t are to depend on its past, a natural generalisation involves replacing the 
constants b and a by some linear functionals on an appropriate function space, say the 
space of continuous functions on a bounded interval. Then we are led by the Riesz' 
representation theorem to the following stochastic differential equation with delay: 

(1) dX(t) = ( f X(t + u)fi(du) )dt+{ I X(t + u)v{du) j dW(t) 

\J[-a,0] J \J[-a,0] J 

for all t ^ and some measures fi, v. We call the solution X of this stochastic differ- 
ential equation geometric Brownian motion with delay and its asymptotic behavior in 
mean square will be characterised in this article. In contrast to the geometric Brow- 
nian motion without delay no explicit representation of X is known from which the 
asymptotic behaviour can be inferred directly. 

Equation (1) is a stochastic functional differential equation with diffusion coeffi- 
cient depending on the past. Such equations may exhibit most irregular asymptotic 
behaviour, see for example Mohammed and Scheutzow [7] for the noisy feedback equa- 
tion. Nevertheless, for a much more general equation than (1) a wide variety of sufficient 
conditions have been established guaranteeing stability in some sense. An exhaustive 
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list of researchers and papers are not quoted here, but a good selection of such results 
are collated in Mao [6] and Kolmanovskii and Myskhis [5] . Despite this activity over the 
last twenty-five years to the best of our knowledge no sufficient and necessary condi- 
tions are known for stability, even for the linear equation (1). By contrast, in this work 
we are able to find necessary and sufficient conditions which characterise completely 
the asymptotic behavior of the solution in mean square. An interesting by-product 
of this stability characterisation is the observation that a deterministic solution may 
transpire to be a solution of the stochastic equation. This feature cannot arise in linear 
non-delay stochastic equations. 

The proof of our stability characterisation relies on the fact that a non-negative 
functional of the process has expected value which satisfies a deterministic linear re- 
newal equation. The asymptotic behaviour of this functional is characterised by the 
renewal theorem; once this characterisation has been obtained, it is straightforward to 
characterise the asymptotic behaviour of the mean square. 



2. Preliminaries 

We first turn our attention to the deterministic delay equation underlying the sto- 
chastic differential equation (1). For a fixed constant a)0we consider the determin- 
istic linear delay differential equation 

x(t) = / x(t + u) /j,(du) for t ^ 0, 

(2) J[-a,0] 

x(t) = ip(t) for t £ [—a, 0], 

for a measure n G M = M[— a, 0], the space of signed Borel measure on [— a,0] with 
the total variation norm IHI^y. The initial function tp is assumed to be in the space 
C[— a, 0] := {V> : [— a, 0] — > 1R : continuous}. A function x : [— a, oo) — > R is called a 
solution of (2) if x is continuous on [— a, oo), its restriction to [0, oo) is continuously 
differentiable, and x satisfies the first and second identity of (2) for all t ^ and 
t € [— a, 0], respectively. It is well known that for every ip G C[— a,0] the problem (2) 
admits a unique solution x = x(-,<p). 

The fundamental solution or resolvent of (2) is the unique locally absolutely contin- 
uous function r : [0, oo) — > 1R which satisfies 

r(t) = l+ [ f r(s + u) dfi(u) ds for t^ 0. 

JO J [max{-«,-s},0] 

It plays a role which is analogous to the fundamental system in linear ordinary differ- 
ential equations and the Green function in partial differential equations. Formally, it is 
the solution of (2) corresponding to the initial function ip = l{o}- For later convenience 
we set r(n) = for u £ [— a,0). 

The solution x(-,<p) of (2) for an arbitrary initial segment ip exists, is unique, and 
can be represented as 



(3) 



x(t, (p) = <p(0)r(t) + / / r(t + s — u) <p(u) du n(ds) for t ^ 0, 
J[-a,Q] Js 
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cf. Diekmann et al [3, Chapter I]. The fundamental solution converges for t — > oo to 
zero if and only if 

(4) vo(jj) := sup \ K(A) : A € C, A- / e Xs fi(ds) = 1 < 0, 

I ■/[-a.o] J 

where 3ft(,z) denotes the real part of a complex number z. In this case the decay is 
exponentially fast (see Diekmann et al [3, Thm. 5.4]) and the zero solution of (2) is 
uniformly asymptotically stable. In this situation we have for every solution x(-,<p) of 
(2): 

\x(t,<p)\ < ce^W for all 0, 

and for a constant c > depending only on the initial function <p. 

Let us introduce equivalent notation for (2). For a function x : [— a,oo) — > R we 
define the segment of x at time t ^ by the function 

Xj : [—a, 0] — ► R, Xt(u) := x(t + u). 

If we equip the space C[—a, 0] of continuous functions with the supremum norm Riesz' 
representation theorem guarantees that every continuous functional F : C[— a,0] — ► R 
is of the form 

Fty) = I ^{u)ii{du) 

J[-a,0] 

for a measure ji G M. Hence, we will write (2) in the form 

x(t) = F(xt) for t ^ 0, xq = ip 

and assume F to be a continuous and linear functional on C[—a, 0]. 

Let us fix a complete probability space (fi, J 7 , P) with a filtration (Ftjt^o satisfying 
the usual conditions and let (W(t) : t > 0) be a Wiener process on this space. We 
study the following stochastic differential equation with time delay: 

dX(t) = F(X t ) dt + G{X t ) dW(t) for t > 0, 
X(u) = (p(u) for u G [—ex, 0], 

where -F and G are continuous and linear functionals on C[— a, 0] for a constant a ^ 0. 
As before, we can write the functional G in the form 

G(V>) = / 4>(u) v{du) for all V G C[-/3, 0] 

J[-a,0] 

and for a measure v G M. We note that assuming the same domain [—a, 0] for the 
arguments of the functionals F and G does not involve any restriction or loss of gener- 
ality. 

For every ip G C[— a,0] there exists a unique, adapted strong solution (X(t,<p) : t ^ 
—a) with finite second moments of (5) (cf., e.g., Mao [6]). The dependence of the 
solutions on the initial condition ip is neglected in our notation in what follows; that 
is, we will write x(t) = x(t,tp) and X(t) = X(t,ip) for the solutions of (2) and (5) 
respectively. 
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By Reifi et al [8, Lemma 6.1] the solution (X(t) : t ^ —a) of (5) obeys a variation 
of constants formula 



(6) 



X(t) 



x(t) + f* r(t - s)G(X s ) dW(s), t ^ 0, 



<f(u), 



u € [-a,0], 



where r is the fundamental solution of (2). It is to be noted that this equation does 
not supply an explicit form of the solution. 



3. Stability 

The asymptotic behavior of the solution X relies on the stochastic convolution inte- 
gral arising in the variation of constants formula (6). Let us define 



(7) 



Y(t) := G{X t ) for t ^ 0, 



such that the stochastic convolution integral is the convolution of the stochastic process 
Y = (Y(t) : t ^ 0) and the fundamental solution r. The following result shows that 
the functional E[Y 2 ] satisfies a linear convolution integral equation. 

Theorem 3.1. Let (X(t) : t ^ —a) be the solution of (5). Then we have for all t ^ 

(8) E \X{t)\ 2 = \x{t)\ 2 + f r\t -s)E \Y{s)\ 2 ds, 

Jo 

where Y , defined by (7), obeys for all t ^ 0: 

(9) E\Y(t)\ 2 = G 2 (x t )+ f G 2 {r t - s )E\Y{s)\ 2 ds. 

Jo 

Proof. The variation of constants formula (6) and Ito's isometry imply the first asser- 
tion. 

Using again the variation of constants formula we obtain by Pubinis's theorem for 
stochastic integrals for t € [0, a]: 

E\Y(t)\ 2 = E\G(X t )\ 2 



= E 



\ X t (u)v(du) + / X t (u)u(du 
J[-a,-t) J[-t,0] 



= E 



-t) 

(p(t + u) v(du) + I x(t + u) v{du) 

-a-t) J[-t,0] 

+ f ( f + r{t + u- s)Y(s)dW(s)^ v{du) 
J[-t,o] \Jo J 

G(x t )+ f If r(t + u- s)u(du) ) Y(s)dW(s) 

Jo \J[s-t,0] J 

= \G(x t )\ 2 + f G 2 {r t ^ s )E\Y( S )\ 2 ds, 
Jo 



= E 
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where we used in the last line r(u) = for u < 0. Setting v([a, b)) = for all a ^ b ^ —a 
enables us to enlarge the integration domain for G such that we can write also for t ^ a: 

2 



E\Y(t)\ 2 = E 



G{x t )+ I ( f + r(t + u- s)Y(s)dW(s)\ v{du) 
J[-t,o] \Jo J 



:,0] 

We can proceed as above to verify that equation (9) is also satisfied for t ^ a. □ 

We next turn to stating and proving our first stability result. In it, the hypothesis 
r € L 2 (R + ) is employed. We remark that this assumption is necessary if E \X 2 {t)\ — > 
as t — > oo. To see this, note by (8) that E \X 2 (t)\ cannot tend to zero as t — > oo if x(t) 
does not tend to zero. But the latter cannot occur if r is not in L 2 (R + ). 

The function s i— > G(r s ), which we denote by G(r»), is square integrable if r € L 2 (IR+) 
and its norm in L 2 (R + ) is given by 

/ roo \ 1/2 

l|G(r.)|| L2(R+) = \J o (G(r s )) 2 dsj . 

This quantity allows to characterise the asymptotic behaviour of the solution for (5): 

Theorem 3.2. // the fundamental solution r is in L 2 (R + ) then the solution (X(t) : 
t ^ —a) of (5) obeys the following trichotomy: 

(a) if \\G(r,)\\ L 2^ + ^ < 1, then there exists k > such 

lim e Kt E\X(t)\ 2 = 0. 

(b) */||G(r.)|| L2(R+) = l, then 

lim E ixmi 2 - ^ G2fe)& ) (/o°° r2 ( 8 ) & ) < 

" JZ°8G*(r a )d8 <0 °- 

(c) i/ ||G r (r,)|| L 2( ]R+ ) > 1, i/ien £/tere exists k > stic/i that 

Proof. Let us introduce the following functions and measures for t ^ 0: 

y(t):=£|y(t)| 2 /(t):=G 2 (x t ) 
5 (t):=G 2 (r t ) ((dt):=g(t)dt. 
Then we can rewrite (9) as the renewal equation 



(10) y(t) 



= /(*) + / - s) C(ds) for all t > 
Jo 



and the three cases (a) to (c) correspond to whether the renewal equation (10) is 
defective, proper or excessive. 

To give the main idea of the proof we establish (a) first without the convergence 
rate. In case (a) the renewal Theorem [1, Thm 3.1.4] implies 

lim y(t) = -JM = 0, 

t-oo 1 - ||G(r.)|| L2(R+) 
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as /(oo) := lim^oo f(t) = due to x(t) — > for t — > oo. Using the notation for y 
introduced above, equation (8) reads 

(11) E\X(t)\ 2 = x 2 (t) + f r 2 (s)y(t- s)ds. 

Jo 

Consequently, as x(t) — > for t — > oo and r G L 2 (R + ) we arrive at 

/•oo 

lim E\X(t)\ 2 = lim / r 2 (s) 1 [0 1 , (s)y(t - s) ds = 

by dominated convergence. 

We now turn to prove the exponential decay. Because ||G r (r,)|| £/ 2( ]R+ ) < 1 there exists 
d > such that 

oo 

e* s ((ds) = 1. 

Moreover, r G L 2 (R + ) implies that there exists 7 > such that r(t) = C(exp(— jt)) 
and x(t) = C(exp(— jt)) as t — > 00. Consequently, we have /(t) = (9(exp(— 27*)) 
and g(t) = C(exp(— 27*)) as t — > 00. Therefore, we can infer by standard methods in 
renewal theory that y(t) = o(exp(— Aci)) for all k < (27 A which leads to 

r-t 

r 2 (s)y{t — s) ds = o(exp(— Ktj). 

Consequently, equation (11) yields the assertion for all k < 27 A 1?. 

In case (b), the renewal Theorem [1, Thm 3.1.5] implies under appropriate conditions 
on / that 



poo 

lim y(t) = (m(C)r 1 / f(s)ds 

t^oo J 



with m(Q := J °° sG 2 {r s ) ds. Note that m(Q is finite as r tends to zero exponentially 
fast and it is non-zero because J °° G 2 {r s ) ds = 1. Since the measure £ is absolutely 
continuous with respect to the Lebesgue measure it is sufficient for the application of 
the renewal Theorem that the function / be in (L 1 nL°°)(]R_|_) and f(t) — ► for t — > 00: 
both of these conditions are satisfied here. By proceeding as in case (a) we obtain 

= ^ (jf /(•)*) (r r2(s)ds ) ■ 

In case (c), the renewal equation (10) is excessive. Then there exists a unique k > 
such that J R e~ KS ((ds) = 1. Furthermore, m K (() := se~ KS ((ds) is non-zero (see 

[1, Remark 3.1.8]) and m K (() is finite as r decays exponentially fast. As in case (b) 
it is sufficient for the application of the renewal Theorem that the function f K with 
f K (t) := e~ Kt f(t) for t ^ tends to zero and is in L 1 (R + ). These conditions are 
satisfied as x(-,<p) decays exponentially fast and hence, the renewal Theorem [1, Cor. 
3.1.9] implies 

lim e~ Kt y(t) = (m^))- 1 / e~ KS f(s) ds. 
Finally, from (8) we have 

% E \X(t)f = e~ Kt \x(t)\ 2 + f e- KS r 2 (s)e- R[ - t - s) y(t - s) ds, 

Jo 



e -nt ; 
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and so, because of the exponential decay of r, we conclude 



Urn e-E\X( t )f = (£ «-/(.) *) «-V(.) *) . 



□ 



Alsmeyer [1] contains a treatment of the renewal equation which covers equations 
with measures. Similar results may be found in Feller [4]. 

As a corollary of Theorem 3.2 we obtain an equivalence between the asymptotic 
behavior of the mean square of the solution X and a condition on the fundamental 
solution r and the diffusion coefficient G. Naturally, this requires that the solution X 
does not reduce to the solution of the deterministic equation (2), for in this case X 
would not provide any information on the diffusion coefficient. We argue below that 
this situation may occur, and must be excluded in the next corollary. This will also be 
illustrated presently in an example. 

Corollary 3.3. Let the fundamental solution r be in L 2 (R + ) and assume that no 
version of the solution X = X(-,<p) of (5) coincides with the deterministic solution 
x = x(-, ip) of (2). Then we have the following: 



lim E\X{t)\ 2 = { 

t— +oo 



l|G(r.)|| £a(]RH . ) <l ) 
c>0 ^ ||G(r.)|| L2(R+) = l, 
oo ||G(r.)|| L2(R+) > 1. 



Proof. We have to show that the constants in Theorem 3.2 in (b) and (c) describing 
the limiting behavior of X are non-zero, which is equivalent to 

/•oo 

(12) / G 2 (x s (-,p))ds^0. 

Jo 

Because s \— > x s (-,ip) and G are continuous operators, equation (12) does not hold if 
and only if G(xt (•, <p)) = for alii ^ 0. In this case x(-,(p) would be also a solution of 
the stochastic equation (5) which contradicts our assumption because of the uniqueness 
of the solution of (5). □ 

Remark 3.4. If the solution x(-,(p) of the deterministic equation (2) also solves the 
stochastic equation (5) it follows by taking expectation and Ito's isometry that 

(13) G(x t (; <p)) = for all t ^ 0. 

Conversely, condition (13) implies that a version of the solution X(-, ip) of (5) coincides 
with x(-, ip). 

Thus, a sufficient condition for the trichotomy only relying on the initial condition 

p is G(ip) ^ 0. 

A more abstract point of view shows that equation (13) holds true if a generalized 
eigenspace N of the deterministic equation (2) is a subset of the kernel ker(G) of 
the diffusion coefficient G. Then, for every ip G N the segment Xt(-,p) is in N and 
consequently in ker(G), so that (13) is satisfied. For details on eigenspaces and related 
results for equation (2) we refer the reader to Diekmann et al [3] and Hale and Lunel 
[2]. A concrete example is given below. 

We emphasise that the forgoing situation in which a deterministic solution may solve 
a non-trivial linear stochastic differential equation is a very specific feature of stochastic 
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functional differential equations, and cannot occur in linear stochastic ordinary differ- 
ential equations. 

Example 3.5. Let us consider the solution (X(t) : t ^ 0) of the simple equation 

(14) dX{t) = bX(t)dt + {cX(t) + dX{t-a))dW(t) for t > 0, 
with b < 0, c, d <G R and a > 0. For the condition in Theorem 3.2 we calculate 

poo 1 

l|G(r.)lli 2(R+ ) = J G 2 (r s )ds = —(c 2 +d 2 + 2cde ba ). 

By using results on deterministic linear difference equations we obtain for a continuous 
function y : [—a, oo) — > R: 

G(y t ) = for all t ^ cy(i) + dy(t - a) = for all t ^ 

<^ y (t) = y(0)e 7 * for all t ^ -a and cd < 

with 7 := — In — . In the case when b 7^ 7 and cd < 0, we obtain that for every initial 
condition <p that the solution X = X(-, p) obeys: 

(15) lim E\X(t)\ 2 = c 2 + + 2cde fea < -25. 

In the case when 6 < and ca! ^ 0, the equivalence (15) also holds. 

In the non-delay case d = the solution X is the geometric Brownian motion and 
(15) reproduces the well-known and easily calculated fact that E\X(t)\ 2 — > if and 
only if c 2 < —2b. In the pure delay case c = and d 7^ we find that d 2 < —2b is 
necessary and sufficient to guarantee E \X(t)\ 2 — > 0. However, although the condition 
on the noise intensities is the same as for geometric Brownian motion, the rate of decay 
to zero is different. 

If c 7^ and d 7^ then the dependence of the stability region on the drift coefficient 
b is described by b < bo < where 60 is the largest real root of 

c 2 + d 2 + 2cde boa + 2b = 0. 

In particular, we observe that while the stability region for (14) is symmetric in c and 
d, it is not symmetric in the sign of cd. 

We finish by pointing out that equation (14) provides an example of the situation 
already mentioned in Remark 3.4 in which a solution of the deterministic equation (2) 
is also a solution of the stochastic equation (5). To see this take for example c = — e, 
d = 1 and b = 7 for some a > 0. Then, for ip(u) = e 7n , u G [— a,0], the solution 
x(t, (p) = e 7 *, t ^ 0, satisfies G(xt(-, </>))= for all i ^ 0. Thus, x(-, y>) is also a solution 
of the stochastic equation and x(t, p) — > for £ — > 00 even though ||G ! (r # )|| L 2( IR+ ) > 1. 
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